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Abstract 

The one-dimensional overdamped Brownian motion in a symmetric periodic potential modulated 
by external time-reversible noise is analyzed. The calculation of the effective diffusion coefficient 
is reduced to the mean first passage time problem. We derive general equations to calculate the 
effective diffusion coefficient of Brownian particles moving in arbitrary super symmetric potential 
modulated: (i) by external white Gaussian noise and (ii) by Markovian dichotomous noise. For both 
cases the exact expressions for the effective diffusion coefficient are derived. We obtain acceleration 
of diffusion in comparison with the free diffusion case for fast fluctuating potentials with arbitrary 
profile and for sawtooth potential in case (ii). In this case the parameter region where this effect 
can be observed is given. We obtain also a finite net diffusion in the absence of thermal noise. For 
rectangular potential the diffusion slows down in comparison with the case when particles diffuse 
freely, for all parameters of noise and of potential. 
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I. INTRODUCTION 



Recently a considerable amount of analysis has been devoted to investigating transport 
of Brownian particles in spatially periodic stochastic structures, such as Josephson junc- 
tions 1], Brownian motors Q and molecular motors 3]. Specifically there has been great 
interest in studying influences of symmetric forces on transport properties, and in calculating 
the effective diffusion coefficient in the overdamped limit in particular P, U, 0, 0, E | ■ Ana- 
lytical results were obtained in arbitrary fixed periodic potential, tilted periodic potentials, 
symmetric periodic potentials modulated by white Gaussian noise, and in supersymmetric 
potentials lllf. The acceleration of diffusion in comparison with the free dif- 

fusion was obtained in Refs. 4, . At thermal equilibrium there is not net transport 

of Brownian particles, while away from equilibrium, the occurrence of a current (ratchet 
effect) is observed generically. Therefore, the absence rather the presence of net flow of 
particles in spite of a broken spatial symmetry is the very surprising situation away from 
thermal equilibrium, as stated in Refs. . Moreover, the problem of sorting of Brownian 
particles by enhancement of their effective diffusion coefficient has been increasingly inves- 
tigated in the last years both from experimental Q ^| and theoretical point of view 
"% [3]. Specifically the enhancement of diffusion in symmetric potentials was investigated 
in Refs. Q,[l4|. 

Motivated by these studies and by the problem of dopant diffusion acceleration in semi- 
conductors physics 16], we try to understand how nonequilibrium symmetrical correlated 
forces influence thermal systems when potentials are symmetric, and if there are new features 
which characterize the relaxation process in symmetric potentials. This is done by using a 
fluctuating periodic potential satisfying the supersymmetry criterion 5], and with a different 
approach with respect to previous theoretical investigations (see review of P. Reimann in 
Ref. [2J). Using the analogy between a continuous Brownian diffusion at large times and the 
"jump diffusion" model [ljj, [ill, we reduce the calculation of effective diffusion coefficient 
D e ff to the first passage time problem. We consider potentials modulated by external white 
Gaussian noise and by Markovian dichotomous noise. For the first case we derive the exact 
formula of D e ff for arbitrary potential profile. The general equations obtained for randomly 
switching potential are solved for the sawtooth and rectangular periodic potential, and the 
exact expression of -D e // is derived without any assumptions on the intensity of driving white 



2 



Gaussian noise and switchings mean rate of the potential. 



II. FAST FLUCTUATING PERIODIC POTENTIAL 



The effective diffusion coefficient in fast fluctuating sawtooth potential was first investi- 
gated and derived in Ref. In papers 17] we generalized this result to the case of arbitrary 
potential profiles. We consider the following Langevin equation 

dx dU (x) 
dt dx 



v{t) + Z (*) 



(i) 



where x(t) is the displacement in time t, £ (t) and rj (t) are statistically independent Gaussian 
white noises with zero means and intensities 2D and 2D„, respectively. Further we assume 

n 

that the potential U (x) satisfies the supersymmetry criterion [5( 

E-U(x) = u(x-^, (2) 
where L is the spatial period of the potential (see Fig. [I]). Following Ref. Q| and because 




FIG. 1: Periodic potential with supersymmetry. 



we have (x (t)) = 0, we determine the effective diffusion coefficient as the limit 

<* 2 (t)> 



D 



eff 



lim 

t^oo 2t 



(3) 



To calculate the effective diffusion constant we use the "jump diffusion" model 



HQ 



x 



n(0,t) 

(*)= 

i=l 



(4) 



where are random increments of jumps with values ±L and n(0, t) denotes the total number 
of jumps in the time interval (0,t). In the asymptotic limit t — > oo the "fine structure" of 



a diffusion is unimportant, and the random processes x (t) and x(t) become statistically 
equivalent, therefore (x 2 (£)) ~ (x 2 (t)). Because of the supersymmetry of potential U (x) 
the probability density reads P (q) = [5 (q — L) + 5 (q + L)] /2. From Eq. (JH) we arrive at 



where r = (tj) is the mean first-passage time (MFPT) for Brownian particle with initial 
position x = and absorbing boundaries at x — ±L. In fluctuating periodic potentials 
therefore the calculation of -D e // reduces to the MFPT problem. Solving the equation for 
the MFPT of Markovian process x(t) we obtain the exact formula for D e ff 



From Eq. @, Deff > D for an arbitrary potential profile U(x), therefore we have always 
the enhancement of diffusion in comparison with the case U (x) = 0. We emphasize that the 
value of diffusion constant does not depend on the height of potential barriers, as for fixed 
potential |8(, but it depends on its gradient U' (x). The dependencies of effective diffusion 
constant D e ff on the intensity D v of the modulating white noise are plotted in Fig. [2] for 
sawtooth, sinusoidal and piece-wise parabolic potential profiles. 
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(5) 
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FIG. 2: Enhancement of diffusion in fast fluctuating periodic potential. 



III. RANDOMLY SWITCHING PERIODIC POTENTIAL PROFILE 



Now we consider Eq. JIJ where r](t) is a Markovian dichotomous noise, which takes the 
values ±1 with switchings mean rate v. Thus, we investigate the Brownian diffusion in a 
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super symmetric periodic potential flipping between two configurations U (x) and — U (x). 
In the "overturned" configuration the maxima of the potential become the minima and vice 
versa. In accordance with Eq. (J2J) we can rewrite Eq. (£[J) as 



and the non-Markovian process x (t) has Markovian dynamics between flippings. Because 
of supersymmetric potential and time-reversible Markovian dichotomous noise the ratchet 
effect is absent: (x) = 0. All Brownian particles are at the origin at t — and the "jump 
diffusion" model (j3J) and (0) is used. The probability density of random increments is the 
same of previous case and the distribution of waiting times tj reads 

w(t) = W+{t)+ o W - {t \ (8) 

where w + (t) and W- (t) are the first passage time distributions for the configuration of the 
potential with 7/(0) = +1 and 77(0) = — 1 respectively. In accordance with Eq. ([5J). r is the 
semi-sum of the MFPTs r + and r_ corresponding to the probability distributions w + (r) and 
W- (r). The exact equations for the MFPTs of Brownian diffusion in randomly switching 
potentials, derived from the backward Fokker-Planck equation, are 

Dt'[ - U' (x) t j + + u(t--t+) = -1, 

Dt j 1 + U'(x)tL + v(t+-t-) = -1, (9) 

where t + (x) and t_(x) are the MFPTs for initial values 77(0) = +1 and 77(0) = —1 respec- 
tively, with the starting position at the point x. We consider the initial position at x = 
and solve Eqs. Q with the absorbing boundaries conditions t± (±£) = 0. Finally we obtain 
the general equations to calculate the effective diffusion coefficient 

o" - f (x) f f (y) e> (y) dy - = x -ip~ , (10) 



'eff 



1 + T~L { l -j)f( x ) e '^) dx 



D 

where 9 (x) = [r + (x) — r_ (x)] /2. Equations ()10j) and (|11)1 solve formally the problem. 



(11) 
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FIG. 3: Switching sawtooth periodic potential. 
A. Switching sawtooth periodic potential 

In such a case (see Fig. EJ) from Eqs. (jlUj) and after algebraic rearrangements, we 
obtain the following exact result 

D eff 2a 2 (1 + fj.) (A„ + \i + /i 2 cosh 2a) 



£> 2a 2 /i 2 (1 + /i) + 2/x sinh 2 a + 4a//(A (U ) sinh a + 8 (A^ 2 ) sinh 2 (o;/2) 
where = 1 — 3/i + 4/icosha, A M i = 7 — fi + 2a 2 /i 2 , and A m2 = 1 — 6/x + /x 2 . Here 



a = a/ (E/D) 2 + z/L 2 / (2D) and /i = vL 2 D / \2E 2 ) are dimensionless parameters, is the 
potential barrier height. The Eq. (|T2"|) was derived without any assumptions on the intensity 
of white Gaussian noise, the mean rate of switchings and the values of the potential profile 
parameters. We introduce two new dimensionless parameters with a clear physical meaning: 
(3 = E/D, and uo = vL 2 /(2D), which is the ratio between the free diffusion time through 
the distance L and the mean time interval between switchings. The parameters a and \l can 



be expressed in terms of (3 and 00 as a = ^/ (5 2 + u, fi = uj/(3 2 . Let us analyze the limiting 
cases. At very rare flippings (u — > 0) we have a ~ /5, \i — > and Eq. (|T2*Jl gives 



Deff . . 



2 



(13) 



D 4sinh 2 (/3/2) ' 

which coincides with the result obtained for the fixed periodic potential. For very 
fast switchings (u — > 00) the Brownian particles "see" the average potential, i.e. 
[U (x) + {—U (x))] /2 = 0, and we obtain diffusion in the absence of potential. If we put in 
Eq. (H2J) a ~ v 7 ^ t 1 + /3 2 / (2a;)] -> 00 and // = u/(3 2 -> cx), we find 

^-1 + ^- (14) 
6 



The normalized effective diffusion coefficient D e ff/D as a function of the dimensionless mean 
rate of potential switching u, for different values of the dimensionless height of potential 
barriers j3, is shown in Fig. 0] A non-monotonic behavior for all values of (3 is observed. 



D eff /D 




FIG. 4: The normalized effective diffusion coefficient versus the dimensionless switchings mean rate 
of potential u) = uL 2 / (2D) for different values of the dimensionless height of the potential barrier. 
Namely (3 = 3, 7, 9, for the curves a, b, and c respectively. 



D e ff/D > 1 for different values above of u. This value decreases with increasing height of 
the potential barrier. In the limiting case of (3 <C 1, we find from Eq. (fT2*|) 

D eff f3 2 ■ [(l + 2cj)cosh2 v ^7-(4cosh v 4J-3)(l+4 v 4;sinh v ^7-2c;)] 



D 



2uj 2 cosh 2\/uj 



(15) 



and for low barriers we obtain the enhancement of diffusion at relatively fast switchings: 
u) > 9.195. For very high potential barriers (f3 — > oo) and fixed mean rate of switchings u, 
we have a ~ [3 — > oo, \i — > 0, and a 2 fi 



a 



w. As a result, we find from Eq. (JT2J) 
vL 2 



eff 



(16) 



We obtained the interesting result: a diffusion at super-high potential barriers (or at very 
deep potential wells) is due to the switchings of the potential only. According to Eq. (|TT)j) the 
effective diffusion coefficient depends on the mean rate of flippings and the spatial period 
of potential profile only, and does not depend on D. The area of diffusion acceleration, 
obtained by Eq. (|12jh is shown on the plane (/3,u) in Fig. [31 as shaded area. This area lies 
inside the rectangle region defined by f3 > and us > 3.5. 
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FIG. 5: Shaded area is the parameter region on the plane (/?, w) where the diffusion acceleration 
compared with a free diffusion case can be observed. 

B. Switching rectangular periodic potential 

For switching rectangular periodic potential represented in Fig. El the main integro- 
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FIG. 6: Switching rectangular periodic potential. 

differential equation (fTU|) includes delta-functions. To solve this unusual equation we use the 
approximation of the delta function in the form of a rectangular function with small width 
e and height 1/e, and then make the limit e —>■ in the final expression. As a result, from 
Eqs. (fTUJl and (fTTj) we get a very simple formula 

D Pff . tanh 2 (/3/2) 



'eff 

D 



1 



cosh {2^/uJ) 



(17) 



We have slowing down of diffusion for all values of the parameters (3 and u. This is because 
in rectangular periodic potential the Brownian particles can only move by thermal force, 
crossing randomly the potential barriers as in fixed potential. The behavior of the normalized 
effective diffusion coefficient D e ff/D as a function of the dimensionless height of the potential 
barrier f3 for different values of the dimensionless mean rate of switchings u is shown in Fig. [3 
The dependence of D e ff/D versus uj for different values of (3 is shown in Fig. [HI For very rare 
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FIG. 7: The normalized effective diffusion coefficient versus the dimensionless height of potential 
barriers (5 = E / D for different values of the dimensionless switchings mean rate uj = uL 2 /(2D). 
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FIG. 8: The normalized effective diffusion coefficient versus the dimensionless switchings mean 
rate uj = uL 2 /(2D) for different values of the dimensionless height of potential barriers (3 = E/D. 



switchings from Eq. ([17]) we obtain the same result for fixed rectangular periodic potential 



D 



D cosh 2 (,3/2) 



(18) 
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In the case of very fast flippings the effective diffusion coefficient, as for sawtooth potential 
(see Eq. (JHJ)), is practically equal to the free diffusion one 

D r 



'eff 



D 

For relatively low potential barriers we get from Eq. (|17|) 



1 -2e- 2 ^tanh 2 (/5/2). (19) 



2 



D 4 cosh (2y/uJ) ' v ' 

Finally for very high potential barriers D e ff depends on the white noise intensity D 

2D . , 

D e „~ - (21) 

1 + COth y/iO 

IV. CONCLUSIONS 



We studied the overdamped Brownian motion in fluctuating supersymmetric periodic 
potentials. We reduced the problem to the mean first passage time problem and derived the 
general equations to calculate the effective diffusion coefficient -D e //. We obtain the exact 
formula for -D e // in periodic potentials modulated by white Gaussian noise. For switching 
sawtooth periodic potential the exact formula obtained for -D e // is valid for arbitrary inten- 
sity of white Gaussian noise, arbitrary parameters of the external dichotomous noise and 
of potential. We derived the area on the parameter plane (/3, us) where the enhancement of 
diffusion can be observed. We analyzed in detail the limiting cases of very high and very low 
potential barriers, very rare and very fast switchings. A diffusion process is obtained in the 
absence of thermal noise. For switching rectangular periodic potential the diffusion process 
slows down for all values of dimensionless parameters of the potential and the external noise. 
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